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ABSTRACT 


Excessive  slamming  is  the  most  common  cause  of  speed  reduction  for 
frigates  and  destroyers  in  heavy  head  seas.  DREA*s  early  interest  in 
frigate/destroyer  slamming  was  limited  to  the  development  of  computer 
programs  for  ship  seakeeping  performance  prediction.  Published  bottom 
slamming  algorithms  were  adopted  and  refined  to  provide  a  slamming 
prediction  capability  for  DREA  software.  In  follow-on  work  described 
herein,  both  two-dimensional  numerical  simulations  and  three-dimensional 
model  tests  were  performed  by  contractors  to  obtain  better  physical  insight 
and  expand  the  empirical  data  base.  The  two-dimensional  numerical 
simulation  results  required  considerable  smoothing  to  reduce  numerical 
noise,  but  still  yielded  form  factors  in  satisfactory  agreement  with 
two-dimensional  theoretical  and  experimental  results.  The  simulations  also 
gave  unique  insight  into  the  girthwise  development  of  slamming  pressures. 
The  three-dimensional  model  tests  indicated  that  pivoted  drop  tests  in 
waves  may  be  substituted  for  conventional  seakeeping  tests,  but  the  general 
applicability  of  the  results  was  limited  by  difficulties  with  measurement 
of  relative  impact  velocity.  Notwithstanding  the  velocity  measurement 
difficulties,  the  three-dimensional  tests  provided  important  information  on 
longitudinal  and  girthwise  pressure  pulse  velocities.  Further  experimental 
work  will  be  required  to  develop  a  usable  experimental  data  base. 
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RESUME 


Un  tapement  excessif  est  la  cause  la  plus  commune  de  la  reduction 
de  vitesse  des  fregates  et  destroyers  par  fortes  mers  debout.  Les 
premieres  recherches  du  CRDA  sur  le  tapement  des  fregates  et  destroyers  se 
limitaient  a  la  mise  au  point  de  programmes  d'ordinateur  pour  la  prevision 
du  rendement  des  navires  en  termes  de  tenue  a  la  mer.  Des  algorithmes 
publies  de  tapement  ont  ete  adoptes  et  raffines  de  maniere  &  permettre  la 
prevision  du  tapement  au  moyen  du  logiciel  du  CRDA.  Lors  des  travaux 
complement a ires  decrits,  des  simulations  numeriques  bidimens ionne lies  et 
des  essais  t  rid  imens  ionne  Is  sur  modeles  ont  ete  effectues  par  des 
entrepreneurs  afin  d'eclairer  1* aspect  physique  et  d'etendre  la  base  de 
donnees  empiriques.  Un  lissage  considerable  des  resultats  de  simulation 
numerique  bidimensionnelle  etait  necessaire  afin  de  reduire  le  bruit 
numerique,  mais  on  a  tout  de  meme  obtenu  des  facteurs  de  forme  concordant 
d'une  maniere  satisfaisante  avec  les  resultats  bidimensionnels  theoriques 
et  experimentaux.  Les  simulations  ont  de  plus  fourni  un  aperfu  inegale  de 
la  propagation  des  pressions  de  tapement  suivant  le  perimetre.  Les  essais 
t ridimensionnels  sur  modeles  ont  6galement  indique  que  1 1  on  pouvait 
substituer  des  essais  de  chute  avec  pivotement  dans  les  vagues  aux  essais 
classiques  de  tenue  a  la  mer,  mais  que  1  *  applicability  generale  des 
resultats  dtait  limitee  par  des  difficultes  de  mesure  de  la  vitesse 
relative  d* impact.  Malgre  ces  difficultes  de  mesure  de  la  vitesse,  les 
essais  tridimensionnels  ont  fourni  des  renseignement s  importants  sur  les 
vitesses  longitudinale  et  suivant  le  perimetre  de  1*  impulsion  de  pression. 
D'autres  travaux  experimentaux  seront  necessaires  pour  l'obtention  d’une 
base  de  donnees  experimentales  utilisable. 
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NOTATION 


A,  B,  C,  K 

Z . . 

h 

h 

k 

k' 


l 

L 

P 

Ph 

PMAX 

Pt 

P(...) 

S 

t 

T 

/V 

V 


V 

vh 


'on 


V, 


ot 


X 

X 


subscripts  for  pressure  transducer  locations 

assumed  speed  of  sound  in  water 

half-siding 

time,  hours 

form  factor 

corrected  form  factor,  per  equation  12 
grid  dimension  for  numerical  simulation 
length  between  perpendiculars 
impact  pressure 

most  probable  extreme  pressure  in  h  hours 
maximum  slamming  pressure 

most  probable  impact  pressure  in  t  seconds 

probability  of  occurrence 

girth 

time,  seconds 
draft 

threshold  velocity 
vertical  ship/section  velocity 
horizontal  component  of  wave  velocity 
initial  impact  velocity 
velocity  component  normal  to  wave 
velocity  component  tangential  to  wave 
vertical  component  of  wave  velocity 
undisturbed  waterline  corresponding  to  y 
x  corrected  for  surge 


v 


X,  Y,  Z 


hull  coordinate  system 
penetration 


y 

a 

a 

6 

0 

y 

x 

5a 


Z>r\ 


P 

aRM 

aRV 

T 


buttock  angle 

parameter  for  extreme  pressure  probability,  Figure  4 

deadrise  angle 

wave  slope 

surge  ratio 

wave  length 

wave  amplitude 

body-fixed  coordinates,  Figure  8 
density  of  water 
RMS  relative  motion 
RMS  relative  velocity 
pitch  angle 


Note :  The  symbols  h  and  a  both  represent  two  parameters,  but  in  each 

case  the  symbol  is  in  conventional  use,  so  the  duplication  was 
retained.  The  appropriate  definition  is  made  clear  in  the  context 
of  each  application. 
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1. 


INTRODUCTION 


For  warships  of  frigate/destroyer  size,  excessive  slamming  is  the 
“"most  common  reason  for  speed  reduc t ion  in” heavy  sea s ;  thus,  a  principal 
goal  in  hull  form  research  and  development  is  the  improvement  of  slamming 
characteristics.  An  implicit  requirement  for  such  improvement  is  the 
development  of  accurate  slamming  prediction  capabilities. 

DREA's  early  interest  in  slamming  was  linked  to  the  development  of 
computer  programs  for  ship  seakeeping  performance  prediction.  These 
programs,  of  the  common  frequency  domain,  linear  strip-theory  type,  adapted 
existing  methods  to  provide  predictions  of  bottom  slam  pressures  and  total 
slam  force. 

Notwithstanding  development  of  these  prediction  capabilities,  it 
was  recognized  that  there  was  still  a  requirement  for  more  fundamental 
information  on  bottom  slamming,  so  further  work  was  conducted,  both 
analytical  and  experimental.  In  the  first  instance,  the  commercially 
available  computer  code  PISCES  2DELK  was  used  to  simulate  drop  tests  of  2-D 
wedges,  truncated  wedges  and  ship  sections.  In  the  second  approach,  ARCTEC 
Canada  Limited  were  contracted  to  conduct  an  experimental  program. 

The  computer  simulations  provided  time  histories  of  local 
pressures,  total  force  and  section  dynamics,  as  well  as  the  development  of 
the  girthwise  pressure  time  history.  After  smoothing  to  reduce  numerical 
noise,  the  data  were  used  to  derive  form  factors,  which  satisfactorily 
agreed  with  two-dimensional  experimental  results. 

The  experimental  work  was  directed  towards  development  of  a  data 
base  containing  slam-related  quantities  for  warship  hull  forms.  Three 
warship  hull  forms,  with  1 U',  ‘normal*  and  *V'  bow  sections  were  used  in 
vertical  and  pivoted  drops,  in  waves  as  well  as  in  still  water.  Seakeeping 
tests  were  also  carried  out  with  the  V  bow  form.  Measurements  included 
local  peak  pressure,  pressure  rise  time,  pulse  duration  and  frequency,  for 
up  to  24  transducers. 

This  memorandum  describes  the  above  work  and  summarizes  earlier 
slamming  research  at  DREA.  The  two-dimensional  simulations  have  provided 
insight  into  the  fundamental  details  of  the  process,  but  it  is  clear  that 
three-dimensional  model  data  are  essential  to  the  achievement  of  realistic 
ship  slam  pressure  predictions.  Experimental  techniques  must  be  refined 
before  a  useful  data  base  can  be  created. 


2.  DREA'S  BACKGROUND  WORK 

The  catalyst  for  slamming  research  at  DREA  was  the  development  of 
an  in-house  computer  program  for  predicting  ship  motions.  PHHS1  (Pitch 
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and  Heave  in  Head  Seas)  employed  conventional  strip  theory  for  the 
calculation  of  vertical  plane  motions  in  unidirectional  head  seas.  This 
coding  was  subsequently  used  together  with  methodology  proposed  by 
Salvesen,  Tuck  and  Faltinsen2  and  enhanced  by  Schmitke3  to  develop  a  6 
-  degree  o f-  -Freed om  p red ic  t-ion  c apab  i-1  i  ty - - - 

The  utility  of  PHHS  was  greatly  increased  by  the  addition  of 
algorithms  for  added  resistance,  relative  motion  corrections  (wave  profile  * 

and  dynamic  swell-up),  deck  wetness,  slamming  pressures  and  human  tolerance 
to  vertical  motion.  Some  of  these  capabilities  were  unique  among 

contemporary  computer  programs.  Relative  motion  corrections  were  of  J 

particular  importance  in  rectifying  strip  theory’s  typically  poor 
prediction  of  deck  wetness;  for  example,  see  Schmitke*  s  contribution  to  the 
discussion  of  Reference  4. 

2.1  Early  Slamming  Algorithms 

The  first  slamming  calculations  were  based  on  the  statistical 
theory  of  Ochi  and  Motter5  together  with  Chuang's6  experimental  impact 
data.  In  the  conventional  manner,  slamming  was  taken  to  be  dependent  upon 
keel  emergence  and  the  reaching  of  a  threshold  relative  velocity  on 
impact.  The  probability  of  keel  emergence  is: 

P(keel  emergence)  =  P(keel)  -  exp  [-T2  /2a2^]  (1) 

where  T  is  draft  and  is  root  mean  square  relative  motion,  corrected 

for  dynamic  swell-up.  Slam  probability  is  then 

P(slam)  -  exp  [-v2  / 2J2Ry]  P(keel)  (2) 

where  crRy  is  RMS  relative  velocity  and  the  slamming  threshold  velocity, 
v,  is  given  by 

v  =  0.0195  /gL/(0.03  cot  8  +  h/B)/2  (3) 

where  L  is  the  LBP,  8  is  local  deadrise  angle,  h  half-siding  and  B  beam. 

The  probability  of  slamming  is  of  rather  academic  interest 
compared  to  the  consequent  slamming  pressures.  As  per  Ochi  and  Motter, 

PHHS  calculated  the  most  probable  peak  impact  pressure  in  t  seconds: 

Pt  =  pka2RV  In  (taRV/2TraRM  P(keel)}  (4) 

7 

The  section  form  factor,  k,  was  obtained  from  a  statistical  fit  to  Chuang's 
data  as  shown  in  Figure  1: 

y 

k  =  1  +  (l-exp(-5  3))  (tt/2  cot  6  )2  (5) 

The  above  expression  is  a  reasonably  good  model  of  the  experimental  data 
for  6  >  5°. 
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2.2  Improved  S landing  Algorithms 

The  early  slamming  algorithms  were  open  to  question  through  their 
use  of  prism  impact  data,  as  suggested  by  comparative  work  by  Ochi  and 
Bonilla-Norat7  and  Schenzle  et  al8  .  The  reply  to  discussion  of 
Reference  5  dealt  with  the  large  magnitude  of  2-D  section  impact  pressure 
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coefficients,  or  form  factors,  at  some  length.  Order  of  magnitude 
differences  between  drop  and  seakeeping  test  derived  form  factors  were 
often  observed,  especially  at  low  deadrise  angles. 

-  From  analysis  of  model  and  full-scale  seakeeping  data,  Ochi  and 

Motter  chose  a  threshold  velocity  of  12  ft/sec  for  a  520  ft  ship.  Through 
Froude  scaling  this  generalised  to 

v  =  0.0927  (gL)1/2  (6) 

This  expression  may  be  criticised  for  its  independence  of  section  / 

shape.  Schmitke9  proposed  an  alternative  derived  from  Ochi's  work: 

v  -  0.377  (gL/k)1/2  (7) 

which  was  obtained  by  multiplying  equation  (6)  by  ( 16. 56/k)1  ^2  ;  16.56 
is  the  k  value  for  Station  3  of  Ochi's  Mariner  hull-form,  used  in  the 
derivation  of  equation  (6).  In  this  way,  a  threshold  impact  pressure  can 
be  defined: 

P  >_  (p/2)k  v2  =  0. 0711pgL  (8) 

Schmitke  noted  that  equations  (7)  and  (3)  produce  comparable  results  for  * 

V-bowed  warships. 

In  modifying  the  algorithms,  the  very  apparent  conservatism  of 
prism  impact  data  was  recognized  and  the  Ochi-Motter  predictor  based  on 
experiments  with  Mariner  models10  was  adopted  for  cases  with  known 
section  offsets.  A  simplified  option  allowing  the  specification  of 
deadrise  angle  only  was  also  developed  using  the  routines  of  Reference  10 
to  generate  form  factors  for  a  series  of  truncated  wedges  (Figure  2)  where 

5°  £  8  £  20° 

T  =  15  ft 

h  =  0.5  ft. 

For  8  less  than  5  degrees,  k  was  taken  as  30.  Since  at  8  of 
21  degrees  the  Ochi-Motter  and  Stavovy-Chuang1 1  methods  produced  similar 
results,  the  results  of  the  two  methods  were  coupled  together  as  shown  in 
Figure  3.  Further  reference  will  be  made  to  this  combination  later.  ,<li 

Slamming  probability  is  of  little  value  itself,  so  the  Ochi-Motter 
method  was  again  used  to  calculate  the  most  probable  extreme  pressure  in  h 
hours  of  operation: 

p^  =  pk  a2 RVln[3600P(keel)h  a^y/(2n,cr^)  ]  (9) 
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t 


Figure  2  Truncated  Wedge  Approximation 


Figure  3  Form  Factor  For  Truncated  Wedges 
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Figure  4  Slamming  Pressure  Probability  Density  Function 


Figure  4  shows  that  the  probability  of  exceeding  is  too  high 
to  be  an  acceptable  design  value,  so  an  extreme  pressure  is  calculated: 

P^(a)  *  pk  ct2 ^yln[3600P(keel)h  o^y/ ( 2JTa  cjrj^)]  (10) 

where  p^Cot)  is  the  extreme  pressure  whose  probability  of  being  exceeded 
in  h  hours  is  a.  Note  that  both  p^  and  p^(a)  are  independent  of 
V*  Reference  4  regards  0.01  as  a  reasonable  value  for  ot. 


3,  IDENTIFICATION  OF  REQUIRED  WORK 

A  number  of  points  became  apparent  in  evaluating  the  literature, 
the  most  important  being: 

-  that  seakeeping  and  drop  tests  predict  inconsistent  slam 
pressures ; 

-  that  full-scale  and  model-scale  seakeeping  tests  yield 
consistent  slam  pressures;  and, 

-  that  the  available  seakeeping-based  algorithms  are  derived  from 
limited  data,  the  bulk  of  which  is  representative  of  merchant 
ship  rather  than  warship  design  practice. 
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In  view  of  these  deficiencies  ARCTEC  Canada  Ltd.  was  contracted  to 
carry  out  a  comprehensive  series  of  experiments  at  model  scale  including 
three-dimensional  drop  tests  and  seakeeping  tests.  The  contract  called  for 

-  a  state  of  the  art  survey  of  experimental  methods  in  slamming; 

-  an  evaluation  of  current  capabilities  of  experimental  facilities 

-  an  experimental  program  to  acquire  slamming  data  for  warship 
hulls ;  and 

-  development  of  a  preliminary  semi-empirical  prediction  method. 


Figure  5  Pressure  Transducer  Locations  on  Bow  Forms 


The  experimental  program  made  use  of  hull  forms  developed  by  the 
National  Research  Council  of  Canada's  Institute  of  Marine  Dynamics  (IMD) 
for  a  fast  surface  ship  methodical  series.12  Normal,  extreme-V  and 
extreme-U  designs  shown  in  Figure  5  were  selected  for  experiments. 

In  parallel  with  the  experimental  work,  numerical  simulations  were 
performed  to  examine  some  of  the  fundamental  features  of  ship  slamming. 


4.  COMPUTER  SIMULATIONS 

4. 1  Background 

In  order  to  gain  a  greater  insight  into  the  physical  parameters 
affecting  the  slamming  process,  a  series  of  computer  simulation  cases  was 
run.  Computer  code  PISCES  2DELK  (a  trademark  of  Physics  International 
Company,  California)  was  selected  for  the  work.  It  is  available  in  Canada 
on  the  CDC  CYBERNET  system.  Individual  cases  were  set  up  and  run  by  CDC 
CYBERNET.  The  resulting  data  were  transferred  to  DREA  for  reduction  and 
analysis. 
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PISCES  2  DELK  is  a  two-dimensional  finite  difference  code  with 
Eulerian,  Lagrangian  and  shell-structure  capabilities.  It  was  designed  to 
simulate  high  velocity  impacts  and  explosions,  i.e.  over  time  scales 
shorter  by  up  to  several. orders  of  magnitude  than  those  characterising  a 
ship  slam.  The  implications  of  this  to  the  present  work  are  noted  in 
Section  4.3.  The  methods  employed  by  PISCES  2DELK  are  described  in  detail 
in  Reference  13. 

4.2  Test  Cases 


Each  case  run  differed  in  geometrical  parameters  and  in  initial 
impact  velocity.  Three,  two-dimensional  geometrical  forms  were  tested: 
simple  wedges,  wedges  truncated  at  the  keel  to  represent  half-siding,  and 
forward  sections  (here  referred  to  as  frigate  sections)  from  the  'normal* 
model  No.  264.  For  wedges,  the  geometrical  parameter  was  deadrise  angle, 

£,  which  could  take  a  value  of  10,  20,  30,  or  40  degrees.  Sections  3,  4, 
and  5  (0  at  FP  and  20  at  AP)  were  used  from  model  264.  Full  scale  impact 
velocity  was  10,  15,  20  or  30  ft/sec. 

4.3  Numerical  Model 

Simulations  were  run  using  full  scale  dimensions.  The  section 
fell  under  the  influence  of  gravity,  simulating  a  free  drop  test.  Figure  6 
illustrates  the  initial  (time  =  zero)  geometry  for  a  simple  wedge.  The 
water,  and  a  void  region  above  it,  were  discretised  in  an  Eulerian  (fixed) 
grid.  The  void  allowed  upward  deformation  of  the  free  surface.  The 
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Figure  6  Initial  Geometry  for  Numerical  Simulation 
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impacting  section  was  assumed  to  be  rigid,  its  boundary  defined  by  a 
Lagrangian  (moving)  grid.  Standard  material  properties  were  modified  to 
ensure  section  rigidity  and  to  reduce  the  speed  of  sound  in  water,  as  noted 
below. 


The  equations  of  motion  were  solved  explicitly  in  time.  In  order 
to  ensure  stability  of  the  solution,  the  time  step  At  was  restricted  to  a 
fraction  of  Jl/c,  where  l  is  the  minimum  cell  dimension  in  the  Eulerian 
grid  and  c  is  the  assumed  speed  of  sound  in  water.  To  achieve  a  reasonable 
value  for  At,  given  l  ,  the  assumed  speed  of  sound  was  reduced  from  the 
true  value,  with  the  limitation  that  the  local  Mach  number  nowhere  exceed 
0.1,  so  that  the  effects  of  compressibility  on  impact  dynamics  could  be 
neglected. 

4.4  Discretisation  Noise 


The  data  were  generally  quite  noisy,  as  illustrated  by  Figure  7. 
Each  curve  on  the  figure  shows  pressure  as  a  function  of  girth  at  a  time 
indicated  by  the  right-hand  vertical  scale.  For  low  deadrise  cases  the 
pressure  distribution  showed  strong  periodicity,  attributed  to  interaction 
between  the  grid  systems.  The  initial  spikes  on  the  keel,  for  all  cases 
with  half-siding,  were  found  to  be  close  to  the  acoustic  pressure  limit, 
pV0c .  The  pressure  data  could  be  smoothed  by  repeated  local  averaging, 
but  this  destroyed  the  peak  structure  of  the  distribution. 

Total  vertical  force  was  smoothed  by  local  averaging  of  the 
cumulative  impulse  a  number  of  times,  and  then  differentiating  the  result. 


Figure  7  Impact  Pressure  Distributions  for  a  20  Degree 
Truncated  Wedged 
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4.5 


Surge  and  Section  Dynamics 


The  parameters  for  a  simple  wedge  during  impact  are^shown  in 
Figure  8.  Surge  is  conveniently  represented  by  the  ratio  y-x/x.  Since 
the  water  surface  was  free  to  deform  in  the  simulations,  surge  could  be 


Figure  8  Impact  Geometry 


readily  estimated  as  illustrated  by  Figure  9  (the  arrows  in  the  figure 
indicate  magnitude  and  velocity  direction  in  the  field).  The  value  of  y 
observed  was  initially  high  but  rapidly  settled  to  a  near  constant  value. 
These  estimates  of  y  for  the  wedge  cases  are  compared  in  Figure  10  with 
the  theoretical  approximation  of  Geers  et  al14  and  with  two  theoretical 
models  of  Chu  and  Abramson15.  Good  correlation  with  the  flat  plate  model 
can  be  seen,  although  selected  experimental  data  cited  in  Reference  15  are 
consistently  lower  than  the  theoretical  models  at  small  deadrise  angles. 

The  surge  ratio  has  some  significance  for  impact  force  as  shown  by 
Figure  11.  The  curves  identified  as  'R-K  simulation*  were  obtained  by 
integrating  the  equations  of  motion  for  a  wedge  using  a  Runge-Kutta 
method.  Two  values  of  y  were  used  in  these  calculations:  the  simulation 
estimate  for  this  particular  case,  1.69;  and  the  average  of  estimates  for 
all  20  degree  wedge  simulations,  1.77.  The  resulting  5  percent  increment 
in  surge  ratio  produces  a  5  percent  increase  in  vertical  force.  The  third 
curve  is  the  result  of  a  PISCES  simulation  with  50  smoothing  iterations 
applied.  Even  this  amount  of  smoothing  does  not  eliminate  discretisation 
noise  entirely. 
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Figure  9  Frigate  Section  3  after  0.367  Seconds: 
Initial  Velocity  10  FPS 


0  10  20  30  40  50 


/SI*) 

Figure  10  Surge  Ratio  for  Simple  and  Truncated  Wedges 
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Figure  11  Time  Histories  of  Vertical  Force  for  a  20  Degree 
Simple  Wedge  with  Initial  Velocity  10  FPS 


4.6  Impact  Pressure  Distribution 

Although,  as  previously  noted,  the  discretisation  employed  in 
these  simulations  was  too  coarse  to  adequately  distinguish  sharp  peaks,  the 
pressure  signatures  and  the  maxima  occurring  after  the  initial  impact  were 
in  general  fairly  well  represented,  even  after  a  number  of  smoothing 
iterations.  To  better  determine  the  initial  peaks,  finer  grids  and  shorter 
time  steps  could  be  used,  but  at  much  greater  computational  cost. 

A  typical  set  of  pressure  distributions  for  a  truncated  wedge  was 
presented  in  Figure  7.  Very  little  difference  between  simple  and  truncated 
wedges  was  observed.  Some  difference  was  seen  in  the  initial  peaks  but 
these  were  heavily  contaminated  with  noise;  after  a  few  time  steps  the 
distributions  were  virtually  indistinguishable.  After  initial  impact,  most 
wedge  cases  approximated  the  classic  Wagner16  pressure  distribution, 
being  fairly  uniform  in  the  region  of  the  keel,  rising  to  a  peak  (the 
amplitude  of  which  depended  on  smoothing)  at  the  waterline. 

Pressure  distributions  for  the  frigate  sections  showed  different 
characteristics  as  typified  by  Figure  12.  Although  initially  like  the 
wedge  pressure  distributions,  the  keel  pressure  falls  more  slowly  with  time 
while  the  waterline  peak  rapidly  disappears.  This  results  in  the  developed 
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pressure  distribution  having  a  triangular  appearance  as  clearly  shown  in 
the  figure.  The  wedges  with  highest  deadrise,  40  degrees,  exhibited 
intermediate  behaviour,  with  the  change  from  a  Wagner  pressure  distribution 
to  triangular  be ing- cons iderab ly~delayed~compared  with- the  frigate  sections-. 


GIRTH  (FT) 

Figure  12  Impact  Pressure  Distributions  for  Frigate  Section  4 


4. 7  Form  Factor 

Because  of  discretisation  noise,  the  pressure  distribution  alone 
would  not  give  a  good  estimate  of  form  factor,  k,  in  the  basic  impact 
equation: 


PMAX  -  <0/2>kV*  (U) 

A  number  of  alternative  methods  to  evaluate  form  factor  were 
tried;  none  were  entirely  satisfactory  for  frigate  sections.  One  of  the 
better  methods  for  wedges  consisted  of  obtaining  a  linear  least-squares 
unbiased  estimate  of  k  from  the  total  vertical  force  produced  by  a  modified 
Wagner  pressure  distribution.  If  k  is  constant,  not  only  should  the 
estimate  be  asymptotic  in  time,  but  also  uncertainty  is  progressively 
reduced.  This  was  found  to  be  so  for  most  wedge  cases,  as  illustrated  by 


13 


Figure  13.  As  shown  in  the  figure,  after  some  initial  excursions,  the 
estimate  stayed  at  a  plateau  for  a  while  before  falling  slowly  in  value. 
This  fall  coincides  with  the  decay  of  the  Wagner  pressure  distribution. 

Form  factor,  taken  to  be  the  aforementioned  plateau  value,  is 
plotted  versus  wedge  deadrise  in  Figure  14.  The  curves  plotted  for 
comparison  are:  Wagner  -  the  classical  result,  Drop  Tests  -  equation  (5), 
and  SC-OM  -  Figure  3. 


Figure  13  Convergence  of  Form  Factor  Estimate  for  Simple 

20  Degree  Simple  Wedge  with  Initial  Velocity  10  FPS 
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Figure  14  Form  Factor  for  Simple  and  Truncated  Wedges 
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5. 


EXPERIMENTAL  WORK 


The  experimental  program  was  preceeded  by  an  extensive  review  of 
the  literature*  In  all  published  sources,  there  is  serious  discrepancy 
between  the  results  of  the  various  types  of  slamming  investigations. 
Reconciling  three-dimensional  drop  and  seakeeping  test  results  is  at  best 
difficult,  so  an  experimental  program  was  derived  that  would  show 
comparison  of  the  results  of  various  techniques  for  the  three  hull  forms 
given  in  Section  3. 

As  planned,  the  test  program's  main  objectives  were: 

-  to  determine  the  form  factor  (or  k-values)  at  each  instrumented 
forward  station; 

-  to  compare  the  measured  form  factors  with  those  predicted  by  the 
Ochi-Motter  (OM)  and  the  Stavovy-Chuang  (SC)  methods; 

-  to  examine  the  effects  of  hull  geometry  on  slam  pressures; 

-  to  develop  a  method  for  calculating  the  spatial  and  temporal 
variation  of  bow  slam  force. 

The  significant  parameters  measured  during  the  experiments  were 
the  peak  pressure  value,  rise  time  (from  zero  to  maximum),  duration  of 
pressure  pulse,  the  time  between  pulses,  and  the  relative  velocity  between 
the  hull  and  the  water  surface  at  the  point  of  contact. 

Three  basic  drop  test  configurations  were  selected,  based  on 
Schenzles'  work8: 

-  even  keel  free  drops  in  calm  water, 

-  pivoted  (rotational)  drops  in  calm  water,  and 

-  pivoted  (rotational)  drops  in  regular  waves; 

as  shown  in  Figure  15.  Additionally ,  seakeeping  experiments  in  regular 
head  waves  were  included  in  the  program. 

5.1  Drop  Test  Programme 

The  drop  tests  were  carried  out  in  a  60  ft  long  basin  with  a 
plunger-type  wave  maker  at  one  end  and  a  sloping  gravel  beach  on  the 
other.  A  modified  towing  carriage  capable  of  accommodating  the  drop  test 
rig  and  achieving  high  forward  speed  was  used  for  the  tests.  The  basin  had 
a  10  ft  wide  wave  channel  in  its  center  with  transparent  walls  at  the  test 
section. 
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( o )  FREE  (EVEN- KEEL )  DROPS  IN  STILL  WATER 
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DROP  HEIGHT  h 


STILL  WATER 


SURFACE 

( b )  PIVOTED  (ROTATIONAL)  DROPS  IN  STILL  WATER 
U 


DROP  ANGLE »- 


(c)  PIVOTED  (ROTATIONAL)  DROPS  IN  REGULAR  WAVES 


Figure  15  Drop  Test  Types 


The  three  models  shown  in  Figure  5  were  fitted  with  twenty-six 
Kulite  XTMS  1-190  series  diaphragm-type  pressure  transducers  effective  in 
the  range  0  to  50  psi.  The  gauges  were  installed  along  the  keel  on  the 
centreline  and  on  one  side  of  the  lower  hull  from  the  F.P.  to  station  7. 
Wherever  possible  three  transducers  were  installed  at  positions  less  than 
1/10  draft  (0 . IT)  from  the  keel,  where  maximum  pressures  were  expected. 
Above  the  0.1T  line,  transducer  locations  were  chosen  to  provide  a 
reasonably  accurate  girthwise  pressure  profile,  and  a  few  were  located  near 
the  load  water  line  forward  to  measure  flare  slamming  effects.  The  model 
was  also  fitted  with  a  motions  package  consisting  of  accelerometers  and 
attitude  indicators. 

The  drop  test  program  was  designed  to  achieve  a  wide  range  of 
relative  velocities  within  an  acceptable  time  and  cost.  A  total  of  849 
drop  tests  were  performed  with  the  three  models.  The  vertical  drop 
velocity  was  readily  controlled  by  the  drop  height,  which  was  restricted  to 
a  maximum  of  12  inches  to  avoid  complete  submergence  of  the  bow.  For  the 
same  reason,  the  maximum  initial  trim  angle  in  the  pivoted  tests  was  14 
degrees.  Calm  water  drops  were  repeated  three  times.  Drops  in  regular 
waves  had  no  repeat  runs  due  to  the  difficulty  in  repeatedly  dropping  the 
model  at  exactly  the  same  position  on  the  wave.  Table  I  gives  the  ranges 
of  the  drop  test  variables. 
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TABLE  I  Ranges  of  Drop  Test  Variables 


Minimum 

Increment 

Maximum 

Drop  Heights 

2" 

2" 

12" 

Initial  Trim  Angles 

2° 

2° 

14° 

Approximate  Forward  Speeds 

0 

0 . 5m/ s 

3m/ s 

5.2  Seakeeping  Test  Programme 

To  complement  the  drop  tests,  ARCTEC  conducted  a  limited  number  of 
seakeeping  tests  using  the  National  Research  Council  (NRC)  of  Canada* s 
facilities  in  Ottawa.  Although  the  original  plan  included  seakeeping 
experiments  with  all  three  models,  experiments  were  carried  out  only  with 
the  V-forra  model  due  to  data  handling  difficulties. 

The  model  was  self-propelled  in  regular  head  waves  in  NRC's  ship 
model  experiment  tank,  with  the  model-carriage  connections  allowing  freedom 
in  pitch,  heave  and  surge  only.  The  model  was  fitted  with  the  same 
pressure  transducer  array  as  in  the  drop  tests,  together  with 
accelerometers  and  sonic  probes  to  derive  motions.  Relative  motions  were 
measured  by  means  of  capacitance  probes  on  the  model  surface.  The 
capacitance  relative  motion  probes  proved  to  be  unsuitable  for  the 
derivation  of  instantaneous  relative  velocities.  Without  relative  velocity 
at  impact,  form  factors  could  not  be  derived  accurately. 

The  seakeeping  experiments  were  performed  in  three  groups.  In  the 
first  case,  Froude  number  and  wave  frequency  were  varied  to  establish 
slamming  zones.  Table  II.  Then,  wave  slope  was  varied  in  the  severe 
slamming  region  as  shown  in  Table  III.  Finally,  the  bracketed  experiment 
in  Table  II  was  repeated  eleven  times  to  obtain  a  measure  of  the  scatter  in 
the  results. 

The  89  test  runs  were  expected  to  generate  about  46  000  pressure 
peaks,  calling  for  an  unreasonable  analysis  effort.  As  a  result,  the 
analysis  was  restricted  to  seven  transducers,  using  a  PDP-11  computer 
digitising  at  20  000  samples  per  second  in  the  region  of  slams.  Even  with 
analysis  of  only  7  of  24  possible  transducers,  50  hours  of  computer  time 
were  required  for  digitising. 
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( a  )  SHIP-WAVE  VELOCITIES  AT  IMPACT 


SHIP  VELOCITIES  Vn  ,Vv 

WATER  VELOCITIES  NORMAL  8  TANGENTIAL 

WAVE  HEIGHT  =  h 
WAVE  LENGTH  -  X 


WAVE  SURFACE  Von  ,V0« 
HORIZONTAL  .VERTICAL  u  ,  v 


SHEER  PLAN 
a  =  BUTTOCK  ANGLE 


(  b)  SHIP  CO-ORDINATES  AND  ANGLES 


Figure  16  Slam  Geometry  Definition 
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PRESSURE (PS  I ) 


5.3 


Results 


5.3.1  Drop  Tests 

For  the~  drop  tests,  the  c  1  ear-wa  Tied ~wave~c"  hartne  1  all  owed  use  o f ~sr 
video  camera  to  record  impacts,  so  that  impact  velocity  and  form  factor 
could  be  reliably  determined.  The  keel  slam  pressure  is  defined  as  follows, 

p  =  (p/2)k'V2  (12) 

k,=  k/(l  +  tan  0  cos20i)  (13) 

V  *  Vh  (Vot  cos  0  +  Von  sin  0)  sin  (a  +  t)  + 

Vv  (Von  cos  ©  -  Vot  s in  0 )  cos  fa  -  T )  (14) 

where  the  angles  are  defined  in  Figure  16  and  the  nomenclature. 


- 95%  CONFIDENCE  LINE 


Figure  17  Typical  Pressure  Plots  from  Stillwater  Pivoted  Drops 
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Figure  18  shows  a  comparison  of  form  factors  derived  from  the  drop 
tests,  together  with  the  SC-OM  line  first  shown  in  Figure  3.  As 
anticipated,  there  are  differences  between  the  flat  and  pivoted  drop  form 
factors,  the  former  giving  clearly  larger  pressures,  except  at  low  deadrise 
angles^ — The~heavy  rat ter- and- reduced  pressures— seen— at— deadrise— angles 
below  30  degrees  could  be  due  to  air  entrapment,  but  it  is  hard  to  credit 
this  at  deadrise  angles  above  five  degrees*  The  pivoted  drop  results  run 
into  the  upper,  Ochi-Motter ,  portion  of  the  SC-OM  line  in  Figure  18,  but 
lie  nearly  an  order  of  magnitude  higher  than  the  lower,  Stavovy-Chuang 
portion.  A  plausable  explanation  for  this  disconcerting  behaviour  is 
available. 

Schmitke9  produced  the  composite  SC— OM  line  by  joining  the  two 
predictions  where  they  crossed.  Unfortunately,  whereas  Ochi— Motter  predict 
k,  Stavovy-Chuang  predict  kj/1*17,  where 

k  -  k^/cos4  8  (15) 

The  chain-dotted  line  in  Figure  18  shows  that  correction  of  the  lower 
portion  of  the  SC-OM  by  the  cos4  8  divisor  would  improve  agreement 
considerably. 

Figure  18  also  shows  that  keel  form  factor  results  for  normal,  U- 
and  V-form  models  collapse  to  a  single  line  when  plotted  to  a  base  of  8  > 
implying  that  hull  form,  per  se,  is  not  so  important  as  deadrise  angle* 

This  is  the  conventional  assumption. 


Figure  19  Comparison  of  Form  Factor  from  Different 
Drop  Test  Types 

Figure  19  compares  derived  form  factors  at  individual  stations  to 
indicate  differences  between  test  techniques.  Even  in  the  reduced  pressure 
zone  (see  Figure  18,  8<25°)  pivoted  drops  in  waves  and  in  calm  water 
yield  consistent  form  factors.  Conversely,  comparison  of  the  pivoted  calm 
and  flat  calm  results  does  not  even  suggest  a  relationship  between  the 
results  of  the  two  techniques.  Perhaps  this  should  be  an  expected  result 
given  the  level  of  scatter  in  the  flat  calm  line  (Figure  18)  at  all 
deadrise  angles.  The  value  of  flat  drops  in  calm  water  must  be  questioned. 
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Velocity  of  Pressure  Pulse  Along  Keel  Model  No.  266  (Extreme  U) 

PIVOTED  DROPS  IN  WAVES 


Model  No.  264  (Normal) 


PIVOTED  DROPS  IN  WAVES 


Model  No.  267  (Extreme  V) 


PIVOTED  DROPS  IN  WAVES 
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40  60 

VELOCITY  (FT/s> 


Figure  20  Slam  Pressure  Pulse  Velocities  cont'd 
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PEAK  SLAM  PRESSURE  AT  KEEL  (PSI ) 


Figure  20  shows  histograms  of  slam  pressure  pulse  velocity  along 
both  the  keel  and  the  girth  for  pivoted  drops  in  waves.  The  average 
pressure  pulse  velocity  is  consistently  higher  than  predicted  by  Ochi  and 
Motter. 

5.3.2  Seakeeping  Tests 

Failure  of  the  seakeeping  model  relative  motion  gauges  prevented 
the  derivation  of  instantaneous  relative  velocities  at  impact,  so  it  was 
necessary  to  use  pitch  and  heave  amplitudes  and  their  average  phases 
relative  to  the  waves  to  estimate  impact  velocities.  This  was  the  source 
of  considerable  scatter  in  pressure  velocity  plots,  so  form  factors  were 
not  derived;  however,  the  seakeeping  data  should  not  be  discounted. 

It  is  revealing  to  superimpose  seakeeping  and  pivoted  drop 
results,  as  in  Figure  21,  for  consistency  of  the  two  methods  is  clear; 
however,  there  is  a  disturbing  amount  of  scatter  at  low  pressures.  There 
is  also  a  suggestion  that  a  common  envelope  encloses  the  two  results;  if 
so,  the  seakeeping  tests  might  have  benefitted  from  more  extreme 
waveheight,  notwithstanding  difficulties  with  relative  velocity 
measurement.  More  important,  a  statistical  estimate  of  the  upper  envelope 
is  potentially  as  useful  a  measure  of  pressure  as  the  mean  line  used 
conventionally. 

Figure  22  shows  the  poor  repeatability  of  peak  slamming  pressures, 
with  near  order  of  magnitude  variations.  The  variations,  probably  due  to 
the  typical  characteristics  of  NRC's  pneumatic  wavemaker,  would  not  be  of 
importance  had  relative  velocity  been  measured,  but  the  absence  of  such 
measurements  prevents  the  derivation  of  reliable  form  factors. 


Figure  21  V-Form  Model  Pressure-Velocity  Relationships : 

Comparison  of  Pivoted  Drop  and  Seakeeping  Tests 
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9  !0  !3  !7  !6  15  !8  !5  17  !4  !!  #  SLAMS 
Figure  22  Repeatability  of  Slam  Pressures  During  Seakeeping  Tests 


6.  S IMPLATION /MODEL  TEST  COMPARISONS 

Direct  comparison  between  simulation  and  model  test  results  is 
precluded  by  the  noise  content  of  the  simulations  and  by  different  impact 
dynamics  in  each  case.  Nonetheless,  this  section  highlights  qualitative 
similarities  observed  in  the  development  of  the  pressure  signature  around 
the  girth. 

A  typical  set  of  model  test  pressure  traces  is  sketched  in  Figure 
23.  These  data  are  at  station  5  for  a  pivoted  drop  in  calm  water  with  an 
initial  full  scale  impact  velocity  of  approximately  20  fps.  Pressures  were 
recorded  on  the  keel  and  at  three  locations  A,B  and  C  towards  the  waterline 
as  shown  in  Figure  5.  Different  levels  of  ambient  noise  on  each  trace 
reflect  different  calibrations  for  each  transducer. 

The  features  seen  in  this  figure  are  typical.  The  large 
oscillations  on  the  keel  trace  may  be  due  to  air  entrapment.  At  locations 
close  to  the  keel,  A  and  B,  the  pressure  rise  is  almost  instantaneous;  at 
location  C  it  is  generally  one  to  several  orders  of  magnitude  slower, 
depending  upon  the  initial  impact  velocity.  It  is  not  possible,  given  the 
•i  spatial  resolution  of  these  pressure  measurements,  to  determine  whether 

this  pressure  rise  character  change  is  gradual  or  abrupt.  Also  note  that 
the  initial  peak  remains  reasonably  constant  between  the  keel  and  B,  but  is 
significantly  reduced  at  C. 

The  dynamics  of  impact  are  different  between  a  pivoted  three- 
dimensional  drop  and  a  simulated  two-dimensional  drop.  Host  notably,  the 
times  at  which  off-keel  peaks  are  recorded  in  the  example  shown  in  Figure 
23  suggest  that  the  model  is  still  accelerating  after  initial  impact, 
whereas  the  equivalent  two-dimensional  simulation  demonstrates  a  small 
deceleration  almost  immediately. 
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Despite  the  above,  and  ignoring  the  keel  oscillations,  the 
pressure  traces  in  Figure  23  are  consistent  with  the  observations  in 
Section  4.6  that  the  girthwise  pressure  distribution  is  initially 
Wagner-like,  that  is,  peaked  at  the  waterline,  (at  least  to  B).  Further 
around  the  girth,  the  pressure  pulse  becomes  much  more  spread  out  resulting 
in  a  more  uniform  pressure  distribution,  tending  to  triangular  as  the  pulse 
disappears  altogether.  In  Figure  24,  data  from  drop  438  are  plotted  with 
data  from  the  corresponding  simulation.  The  development  of  the  pressure 
distribution  as  described  above  can  be  seen  in  both  cases. 

The  transition  from  a  Wagner  to  a  triangular  pressure  distribution 
occurs  as  local  deadrise  is  increasing,  and  may  be  associated  with  a 
critical  value  of  deadrise.  The  simulations  suggest  that  this  transition 
is  fairly  abrupt  in  time,  as  shown  in  Figure  24  between  step  72  and  step 
84.  The  value  of  deadrise  for  this  transition,  estimated  from  the 
simulation  data,  is  shown  in  Figure  25.  There  is  no  systematic  variation 
of  transition  position  with  impact  velocity. 

The  width  of  the  transition  region  in  Figure  25  is  due  to 
numerical  scatter  in  the  simulations.  Since  the  critical  deadrise  angle 
for  transition  differs  for  each  station  in  Figure  25,  transition  is  not  a 
function  of  deadrise  angle  alone,  unlike  form  factor.  In  addition, 
extrapolation  suggests  that  for  fuller  sections,  the  critical  value  of 
deadrise  may  be  quite  low. 


40 

20 

co  0 

o. 

UJ 

20 

CO 

CO  - 

uj  0 

<r 

Q. 

<  20 

-I 

40  0 
0 

0  50  100  150  200 

MILLISECONDS  FROM  KEEL  IMPACT 

Figure  23  Drop  438, Sketch  of  Pressure  Signatures  at 
Station  5,  to  full  scale 


28 


Figure  24  Comparison  of  Pressure  Distribution  Development 
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Figure  25  Transition  from  Wagner  to  Triangular  Girthwise 
Pressure  Distribution  from  Simulation  Data 


7.  CONCLUDING  REMARKS 

The  work  described  in  this  memorandum  has  been  aimed  towards  the 
eventual  development  of  accurate  slamming  prediction  capabilities  for 
warship  hulls*  The  initial  tasks,  presented  here,  consist  of  a 
comprehensive  set  of  model  experiments,  undertaken  by  ARCTEC  Canada  Ltd*  , 
to  provide  a  database  of  slamming  data;  and  a  series  of  numerical 
simulations  to  provide  further  insight  into  the  fundamental  processes  of 
slamming  and  to  aid  with  interpretation  of  the  data*  The  simulations  were 
restricted  by  computational  complexity  to  two  dimensions* 

Results  from  the  numerical  simulations  agree  well  with 
two-dimensional  theories  and  drop  test  data*  The  surge  ratio,  |i »  plays 
a  significant  role  in  two  dimensional  drop  dynamics,  although  it  is  not 
clear  to  what  extent  this  translates  into  three  dimensions.  Numerical 
noise  made  it  necessary  to  smooth  the  results  and  made  peak  pressure 
determination  uncertain. 

A  notable  aspect  of  the  simulation  results  was  the  transition  of 
the  developing  bottom  pressure  distribution'  from  a  classic,  Wagner,  form 
to  a  uniform,  or  triangular,  shape.  This  transition  appears  to  be  fairly 
abrupt,  and  may  be  associated  with  quite  low  local  deadrise  angles  for 
full  sections,  although  it* is  clearly  not  a  function  of  deadrise  alone. 
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The  transition  was  also  noted  on  the  higher  deadrise  angle  wedge 
simulations,  but  was  much  delayed.  Although  spatial  resolution  of  the 
bottom  pressure  was  low,  the  experimental  data  appears  to  confirm  the 
existence  of  this  phenomenon  for  three  dimensional  hulls. 

In  the  experimental  program,  the  results  from  the  extensive  drop 
test  series  and  the  more  limited  seakeeping  experiments  demonstrated  that 
pivoted  drop  testing  is  a  useful  model  experimental  technique  from  which 
the  slamming  characteristics  of  warship  hull  forms  can  be  evaluated.  A 
good  consistency  was  found  between  pivoted  drops,  both  in  calm  water  and 
in  waves,  and  seakeeping  tests. 

The  experimental  results  show  a  fair  degree  of  scatter, 
particularly  at  low  impact  velocities.  These  uncertainties  are  primarily 
in  the  measurement  of  impact  velocity,  which  needs  further  attention  in 
future  work.  Nonetheless,  predictions  of  pressure-velocity  relationships 
and  the  derived  form  factors  are  satisfactory  for  design  applications. 

A  notable  departure  from  previous  theories  is  that  the  slam 
pressure  pulse  velocities  observed  for  pivoted  drops  in  waves  were 
consistently  higher  than  predicted  by  Ochi  and  Motter.  Also,  the 
previously  adopted  Stavovy  Chuang-Ochi  Motter  line  appears  to  underpredict 
form  factors  at  deadrise  angles  greater  than  25  degrees,  apparently  due  to 
misinterpretation  of  the  original  Stavovy-Chuang  data. 

Consideration  of  the  progress  achieved  this  far  leads  to  the 
following  recommendations  for  further  investigation: 

-  An  improved  technique  for  impact  velocity  measurement  is 
required,  possibly  using  some  form  of  emergence  probe. 

-  Specific  experiments  and  numerical  simulations  should  be 
conducted  to  clarify  further  some  of  the  fundamental  features  of 
slamming,  such  as  the  spatial  distribution  and  time  history  of 
pressures,  particularly  in  three  dimensions. 

-  A  larger  experimental  database  is  required  to  allow  the 
establishment  of  reliable  empirical  prediction  methods. 
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